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Abstract 

The spherically symmetric thin shells of the barotropic fluids with the linear equation of 
state are considered within the frameworks of general relativity. We study several aspects 
of the shells as completely relativistic models of stars, first of all the neutron stars and 
white dwarfs, and circumstellar shells. The exact equations of motion of the shells are 
obtained. Also we calculate the parameters of the equilibrium configurations, including the 
radii of static shells. Finally, we study the stability of the equilibrium shells against radial 
perturbations. 
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In present paper we consider some infinitely thin (singular) shells as models of the physical 
entities, which thickness is negligible in comparison with a circumference radius (e.g., the 
circumstellar shells or surfaces of phase transitions). Besides, one may study the shells 
as the simplest (nevertheless useful and instructive) models of "thick" objects, e.g., stars. 
Geometrically the shell is described by a three-dimensional closed hypersurface, embedded 
in the four-dimensional spacetime and dividing it into the two domains, the external (S + ) 
and internal (S~) spacetime. Since the classical works [5], || the theory of surface layers has 
been widely considered in the literature (see ref. Q for details) including modern applications 
and developments in fundamental, quantum and phenomenological directions, see ref. j| and 
references therein. We only point out some questions of principle now. One of them is the 
essential difference between the description of the boundary surface (for instance, surface of 
star) and shell. The boundary surface is the discontinuity of the first kind (mass density 
has a finite jump across the surface) and is described by the Lichnerowicz-Darmois junction 
conditions (the first and second quadratic forms are continuous on the surface). The singular 
thin shell is the discontinuity of the second kind (the density has the delta-like singularity 
on the shell) and is described by the Lichnerowicz-Darmois- Israel junction conditions: the 
first quadratic form (the metric) is continuous, the second one (the extrinsic curvature) has 
a finite jump. 

This circumstance says we must take care at physical interpretation (and application) 
of the shells. For instance, the surface of a star should be described by a boundary surface 
not a singular thin shell. Latter can be used as a model: (a) of the star as a whole system 
||], we replace a thick layer of matter (which has an infinite number of degrees of freedom) 
by the singular shell which has a finite number of collective degrees of freedom (spherically 
symmetric thin shells have only radial degree of freedom), (b) of the stellar internal structure, 
thereby a star can be considered as a set of an infinite number of embedded shells, (c) of the 
circumstellar shells as a set of a finite number of shells. The aim of this paper is to study all 
appropriate cases for the simplest class of equations of state (EOS). 

So, one considers a thin shell with the surface stress-energy tensor of a perfect fluid in 
the general case 

Sab = CTU a U b + p(u a U b + (3) g afe ), (1) 

where a and p are the surface energy density and pressure respectively, u a is the timelike unit 
tangent vector, ^g a b is the 3-metric of a shell's surface. We suppose metrics of the spacetimes 
outside S + and inside X - of a spherically symmetric shell to be of the Schwarzschild form, 
i.e., 

2M±1 'dr' + rW, (2) 



d4 



l _2M± 
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where dft 2 is the metric of the unit 2-sphere. It corresponds to the neutral shell surrounding a 
neutral body of mass M_, and M + appears to be the external observable total mass (energy) 
of the whole system. We assume both masses to be non-negative (though a theory of quantum 
gravity could give rise to spacetimes of negative mass, see ref. W). It is possible to show 
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that if one uses the proper time r then the 3-metric of a shell is 

( 3 )d s 2 = -dr 2 + R 2 dQ 2 , (3) 

where R = R(r) is a proper radius of a shell. Also the energy conservation law for matter 
on a shell (which is an integrability condition for the Einstein equations on a shell) can be 
written as 

d(aWg)+pd(Wg)=0, (4) 



where ^g = J — det (^g a b) = R 2 sin9. In this equation, the first term corresponds to a 
change in the shell's internal energy, the second term corresponds to the work done by the 
shell's internal forces. 

Imposing junction conditions across a shell, the Einstein equations yield the equation of 
shell's motion in the form 



■ 2M+ / •„ 2M_ m 

. 1 + R2 -nr- £ -V 1 + R2 - — = -ir <5 > 

where 

m = AttctR 2 (6) 
is interpreted as an effective rest mass, R = dR/dr is the proper velocity of a shell, e± = 



sign f y 1 + R? — 2M±/Rj . Equations @ - @ together with an equation of state and choice 

of the signs e±, completely determine the motion of the fluid Schwarzschild shell in general 
relativity. Therefore, we must resolve the next two problems, viz., the choice of e± in eq. (||), 
and the choice of EOS. 

(i) The choice of signs e±. It is well-known, that e = +1 if R increases in the outward 
normal direction to the shell ( e.g., it takes place in a flat spacetime), and e = — 1 if R 
decreases (semiclosed world). Thus, only under the additional condition e + = e_ = 1 we 
have an ordinary shell [Q, ||]. In the paper we model ordinary stars not the wormhole-stars, 
therefore, in (||) we suppose e + = e_ = 1. 

(ii) The choice of the equation of state. We confine our discussion to models calculated 
with ultra-high density barotropic state equations of the linear form 

% = ( 7 - 1) (7) 

For physical purposes the cores with 1 < 7 < 2, e.g., 7=2, 4/3, 6/5, 14/13 were considered 
P] as models of neutron stars. They showed a maximum mass of approximately 0.7 M@ 
corresponding to a radius around four Schwarzschild radii. For two-dimensional matter 
eq. (^) implies 

V = W, (8) 

where the reduction of dimensionality of the stress-energy tensor should be taken into ac- 
count, i.e., 2r] = 3(7 — 1). This equation includes the most studied cases: the dust shell p = 
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|, 0, radiation fluid shell a - 2p = [|TTJ, and bubble <r + p = @, g. If rj > it can 
be interpreted as a square component of the vector of a speed of sound in the shell. Then 
for a spatially two-dimensional homogeneous fluid the square speed of sound is 2r/. From the 
physical viewpoint some ry appear to be inadmissible. For instance, if the fluid on a shell is 
required to satisfy the dominant energy condition, a > \p\, one obtains the constraint 

H < 1- (9) 

If a fluid is required to satisfy the causality condition, we get the constraint 

ry < 1/2, (10) 

where one took into account spatial two-dimensionality of the fluid. Nevertheless, the aim of 
this paper is to study the general case of arbitrary 77 [jOJ. So, solving the differential equation 
(Q) with respect to a, we obtain 

a = £LR-*<*+V, (11) 

47T 

where C is the integration constant determined by the specific shell matter. The value of C 
is closely related to the value of surface mass density (or pressure) at fixed R. We consider 
ordinary shells not wormholes, therefore, we require a > 0, hence C > 0. Also it should 
be noted that from (|6|) at positive densities a it follows that M + > M_ for any R and 
R. Otherwise, matching of the spacetimes (^) is impossible. Thus, we have all necessary 
equations (|5|), (^) and (pd|) to study the shells as models of stars. 

First of all we study the static shells. Taking into account (El), the equilibrium conditions 



R = R = |14| for our case read 

k- — k + = 47T<Toro, (12) 

Mr M_ 



471(27? + l)a rl (13) 



where k± = a/1 — 2M±/tq (k + < k-), r$ and cjq are the shell radius R and density at 
equilibrium respectively, o"o is given by (|ll]) at R = 7*0 . Taking into account (|l^), the 
equation ( |i~3P can be rewritten in the nice form 

1 

4?y + 1. (14) 



k+k- 



We can use the system (|l~2j), ( |i~4| ) instead (|i~2"D, (13). The equation (|14|) is very helpful because 
it already gives us both the radius rg and constraint for 7/ independently of the arbitrary 
constant C. 

Proposition. Let we have the neutral ordinary shell of the matter with the EOS (|8|), k± 
are real and positive, M + 7^ M_. Then the shell can have static states only at 

T] > 0. 
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Proof: From eq. (14) it follows that Arj + 1 > 0, and 



(iy + l)~ 2 = 2Af + 
l-2M_/r r 

or 

(4?] + l)- 2 < 1. 

These conditions give the desired inequality rj > 0, Q.E.D. 

Hence it follows that neither the dust shells (rj = 0), nor bubbles (rj = —1) can have static 
states with mass M + . In other words, they should not be considered as static sources of the 
Schwarzschild metric @j. By means of this proposition we can prove the next theorem. 

Theorem. Let the conditions of the previous proposition are valid. Then the radius of 
static shells is given by the expression 



1 

ro = ~ 
Q 



M + + M_ + JAM 2 + 4M+M_(1 - q) 



(15) 



where g = 1 - (4r? + 1)~ 2 , AM = M+ - M_. 

Proof: Perform the conformal transformation £ = ro/2M + , a = M_/M + . The case a = 
corresponds to the hollow shell (M_ = 0), but a = 1 describes the limit case M± — > oo, 
AM — > (infinitely light shells), not the case M + = M_ which for ordinary shells is forbidden 
as trivial (|5|). Therefore < a < 1. Then, squaring eq. (M), we obtain the square equation 
with respect to £ 

gf - (1 + a)f + a = 0, 

which has the two roots 



2£>£± = 1 + a± y(l + a) 2 - 4ag, (16) 
one of them is superfluous, as can be directly verified. The preceding proposition yields 



< g < 1, J (1 + a) 2 - Aag < 1 + a. 
Further, taking into account eq. (|i~6|), the condition £ > 1 can be written as 



± J (I + a) 2 -\ag > 2g- 1 - a. (17) 
If we define the auxiliary function 



z±(o) = 1 + a - 2g ± J (1 + a) 2 - 4ag, 
then we must require z±(g) > 0. It can easily be seen that 

z + (g ^ 1) > 0, z+(l) = 0, z^(g + 0) < 0, z_(0) 
therefore, only the root £+ is admissible, Q.E.D. 
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The radius ro given by eq. ( |l5| ) is illustrated in Fig. |l] for the most physically admissible 
range, (dust) < r\ < 1/2 (ultrarelativistic fluid). For clarity only the two limit cases, a = 
and a = 1, are represented, the curves for other < a < 1 lie between them. The parameter 
1 — a = AM/M+ determines the proper total mass of the shells. As it was mentioned above, 
we can model both the whole star and stellar internal structure. In this connection the model 
with a = (the hollow shell) appears to be a model of the star as a whole system, whereas 
the shells at a / can be models both inhomogeneous stellar structure (the layers with 
smaller r\ lie above those with greater rj) and the circumstellar shells. 

Another important aspect we must consider is the mechanical stability of the shells. Here 
we apply the method developed in ref. [14|]. Redefining the constant C in (pT|), we suppose 



(18) 



CTq 47T 



where z = R/tq, the constants do and vq should be calibrated by means of the system 
(|T2|), (|l4|). It is well-known that the shell's equation of motion can be rewritten in a more 
convenient form. Squaring (||) twice, one obtains 



R 2 + V(R) = 0, 



(19) 



where the potential V(R) is 



V(R) = 1 



2M 4 



R\ 2r > +1 Ak fr \ 2r i +1 

r W ~ T \r) 



where a = AM/Ak, Ak = k + — k-. Calibrating these expressions, we obtain 



r z z 2 + V(z) = 0, 



(20) 



where 



V(z) 



1 + 



g{k\ + l)z 4v+1 + gk\-l 



1 2 



2gk + z 



2r,+l 



The potential satisfies V(z = 1) = dy/dz| 2= i = 0, as can be directly verified. Therefore, 
the stability of the equilibrium is equivalent to V(z) having a minimum at z = 1. By 
differentiation one obtains 



d 2 V 



dz 2 



2=1 



fc|(4r? + l) 2 - 2*4 + 1 
2fci ' 



(21) 



From eq. ( |2l[) and from the proposition proved above it can readily be seen thatd^y/dz 2 | z= i < 
0. Therefore, the equilibrium of the neutral shells with the linear EOS (unlike those with 
the linearized EOS studied by Brady, Louko and Poisson [jl4|) is unstable against radial 
perturbations. In application to models of neutron stars it means that the EOS d8|) can not 
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purely describe them. Apparently only the average layers of the star core can be "made" 
from the matter (||). The central part of the core (the so called "core of the core") and core's 
crust must be consisting of other kind substance to prevent from collapse into black holes 
and explosion respectively. For instance, Witten IE] argued that the strange matter can 
be formed at the quark-hadron phase transition in central part of neutron stars. As for the 
crust that it can be even crystallic |16| . 

Thus, in the present paper some aspects of dynamics of the spherically symmetric shells 
with the linear EOS were considered. Finally it should be noted that the studied shells 
may also describe, besides static (or having adiabatically slow radial motion) circumstellar 
shells and stars, the transient astrophysical phenomena, e.g., the blast waves at supernova 
explosion. Thereby the obtained models of these phenomena will be completely relativistic 
that is essential at high velocities and powerful fields when the spacetime can not be supposed 
flat. 
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Figure 1: Plot of the static shell radius (in units of the Schwarzschild radius) vs rj. The dashed 
curve describes hollow shells (a = 0), the solid curve does the asympotical case a = 1 (infinitely 
light shells). The shaded region corresponds to black holes. 
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